Comprehensive molecular profiling of advanced/metastatic olfactory neuroblastomas.

Olfactory neuroblastoma (ONB) is a rare, locally aggressive, malignant neoplasm originating in the olfactory epithelium in the nasal vault. The recurrence rate of ONB remains high and there are no specific treatment guidelines for recurrent/metastatic ONBs. This study retrospectively evaluated 23 ONB samples profiled at Caris Life Sciences (Phoenix, Arizona) using DNA sequencing (Sanger/NGS [Illumina], n = 15) and gene fusions (Archer FusionPlex, n = 6), whole genome RNA microarray (HumanHT-12 v4 beadChip, Illumina, n = 4), gene copy number assays (chromogenic and fluorescent in situ hybridization), and immunohistochemistry. Mutations were detected in 63% ONBs including TP53, CTNNB1, EGFR, APC, cKIT, cMET, PDGFRA, CDH1, FH, and SMAD4 genes. Twenty-one genes were over-expressed and 19 genes under-expressed by microarray assay. Some of the upregulated genes included CD24, SCG2, and IGFBP-2. None of the cases harbored copy number variations of EGFR, HER2 and cMET genes, and no gene fusions were identified. Multiple protein biomarkers of potential response or resistance to classic chemotherapy drugs were identified, such as low ERCC1 [cisplatin sensitivity in 10/12], high TOPO1 [irinotecan sensitivity in 12/19], high TUBB3 [vincristine resistance in 13/14], and high MRP1 [multidrug resistance in 6/6 cases]. None of the cases (0/10) were positive for PD-L1 in tumor cells. Overexpression of pNTRK was observed in 67% (4/6) of the cases without underlying genetic alterations. Molecular alterations detected in our study (e.g., Wnt and cKIT/PDGFRA pathways) are potentially treatable using novel therapeutic approaches. Identified protein biomarkers of response or resistance to classic chemotherapy could be useful in optimizing existing chemotherapy treatment(s) in ONBs.

The Genetic Architecture of Coronary Artery Disease: Current Knowledge and Future Opportunities.

PURPOSE OF REVIEW: We provide an overview of our current understanding of the genetic architecture of coronary artery disease (CAD) and discuss areas of research that provide excellent opportunities for further exploration. RECENT FINDINGS: Large-scale studies in human populations, coupled with rapid advances in genetic technologies over the last decade, have clearly established the association of common genetic variation with risk of CAD. However, the effect sizes of the susceptibility alleles are for the most part modest and collectively explain only a small fraction of the overall heritability. By comparison, evidence that rare variants make a substantial contribution to risk of CAD has been somewhat disappointing thus far, suggesting that other biological mechanisms have yet to be discovered. Emerging data suggests that novel pathways involved in the development of CAD can be identified through complementary and integrative systems genetics strategies in mice or humans. There is also convincing evidence that gut bacteria play a previously unrecognized role in the development of CAD, particularly through metabolism of certain dietary nutrients that lead to proatherogenic metabolites in the circulation. A major effort is now underway to functionally understand the newly discovered genetic and biological associations for CAD, which could lead to the development of potentially novel therapeutic strategies. Other important areas of investigation for understanding the pathophysiology of CAD, including epistatic interactions between genes or with either sex and environmental factors, have not been studied on a broad scope and represent additional opportunities for future studies.

Multi-platform molecular profiling of a large cohort of glioblastomas reveals potential therapeutic strategies.

Glioblastomas (GBM) are the most aggressive and prevalent form of gliomas with abysmal prognosis and limited treatment options. We analyzed clinically relevant molecular aberrations suggestive of response to therapies in 1035 GBM tumors. Our analysis revealed mutations in 39 genes of 48 tested. IHC revealed expression of PD-L1 in 19% and PD-1 in 46%. MGMT-methylation was seen in 43%, EGFRvIII in 19% and 1p19q co-deletion in 2%. TP53 mutation was associated with concurrent mutations, while IDH1 mutation was associated with MGMT-methylation and TP53 mutation and was mutually exclusive of EGFRvIII mutation. Distinct biomarker profiles were seen in GBM compared with WHO grade III astrocytoma, suggesting different biology and potentially different treatment approaches. Analysis of 17 metachronous paired tumors showed frequent biomarker changes, including MGMT-methylation and EGFR aberrations, indicating the need for a re-biopsy for tumor profiling to direct subsequent therapy. MGMT-methylation, PR and TOPO1 appeared as significant prognostic markers in sub-cohorts of GBM defined by age. The current study represents the largest biomarker study on clinical GBM tumors using multiple technologies to detect gene mutation, amplification, protein expression and promoter methylation. These data will inform planning for future personalized biomarker-based clinical trials and identifying effective treatments based on tumor biomarkers.

Expanding role of gut microbiota in lipid metabolism.

PURPOSE OF REVIEW: This article highlights recent advances in the emerging role that gut microbiota play in modulating metabolic phenotypes, with a particular focus on lipid metabolism. RECENT FINDINGS: Accumulating data from both human and animal studies demonstrate that intestinal microbes can affect host lipid metabolism through multiple direct and indirect biological mechanisms. These include a variety of signaling molecules produced by gut bacteria that have potent effects on hepatic lipid and bile metabolism and on reverse cholesterol transport, energy expenditure, and insulin sensitivity in peripheral tissues. Additionally, host genetic factors can modulate the abundance of bacterial taxa, which can subsequently affect various metabolic phenotypes. Proof of causality for identified microbial associations with host lipid-related phenotypes has been demonstrated in several animal studies, but remains a challenge in humans. Ultimately, selective manipulation of the gut microbial ecosystem for intervention will first require a better understanding of which specific bacteria, or alternatively, which bacterial metabolites, are appropriate targets. SUMMARY: Recent discoveries have broad implications for elucidating bacterially mediated pathophysiological mechanisms that alter lipid metabolism and other related metabolic traits. From a clinical perspective, this newly recognized endocrine organ system can be targeted for therapeutic benefit of dyslipidemia and cardiometabolic diseases.

Multiplatform molecular profiling identifies potentially targetable biomarkers in malignant phyllodes tumors of the breast.

Malignant phyllodes tumor is a rare breast malignancy with sarcomatous overgrowth and with limited effective treatment options for recurrent and metastatic cases. Recent clinical trials indicated a potential for anti-angiogenic, anti-EGFR and immunotherapeutic approaches for patients with sarcomas, which led us to investigate these and other targetable pathways in malignant phyllodes tumor of the breast. Thirty-six malignant phyllodes tumors (including 8 metastatic tumors with two cases having matched primary and metastatic tumors) were profiled using gene sequencing, gene copy number analysis, whole genome expression, and protein expression. Whole genome expression analysis demonstrated consistent over-expression of genes involved in angiogenesis including VEGFA, Angiopoietin-2, VCAM1, PDGFRA, and PTTG1. EGFR protein overexpression was observed in 26/27 (96%) of cases with amplification of the EGFR gene in 8/24 (33%) cases. Two EGFR mutations were identified including EGFRvIII and a presumed pathogenic V774M mutation, respectively. The most common pathogenic mutations included TP53 (50%) and PIK3CA (15%). Cases with matched primary and metastatic tumors harbored identical mutations in both sites (PIK3CA/KRAS and RB1 gene mutations, respectively). Tumor expression of PD-L1 immunoregulatory protein was observed in 3/22 (14%) of cases. Overexpression of molecular biomarkers of increased angiogenesis, EGFR and immune checkpoints provides novel targeted therapy options in malignant phyllodes tumors of the breast.

Mutations in the Kinase Domain of the HER2/ERBB2 Gene Identified in a Wide Variety of Human Cancers.

The HER2 (official name ERBB2) gene encodes a membrane receptor in the epidermal growth factor receptor family amplified and overexpressed in adenocarcinoma. Activating mutations also occur in several cancers. We report mutation analyses of the HER2 kinase domain in 7497 histologically diverse cancers. Forty-five genes, including the kinase domain of HER2 with HER2 IHC and dual in situ hybridization, were analyzed in tumors from 7497 patients with cancer, including 850 breast, 770 colorectal, 910 non-small cell lung, 823 uterine or cervical, 1372 ovarian, and 297 pancreatic cancers, as well as 323 melanomas and 2152 other solid tumors. Sixty-nine HER2 kinase domain mutations were identified in tumors from 68 patients (approximately 1% of all cases, ranging from absent in sarcomas to 4% in urothelial cancers), which included previously published activating mutations and 13 novel mutations. Fourteen cases with coexisting HER2 mutation and amplification and/or overexpression were identified. Fifty-two of 68 patients had additional mutations in other analyzed genes, whereas 16 patients (23%) had HER2 mutations identified as the sole driver mutation. HER2 mutations coexisted with HER2 gene amplification and overexpression and with mutations in other functionally important genes. HER2 mutations were identified as the only driver mutation in a significant proportion of solid cancers. Evaluation of anti-HER2 therapies in nonamplified, HER2-mutated cancers is warranted.

Epigenome-wide association of liver methylation patterns and complex metabolic traits in mice.

Heritable epigenetic factors can contribute to complex disease etiology. Here we examine the contribution of DNA methylation to complex traits that are precursors to heart disease, diabetes, and osteoporosis. We profiled DNA methylation in the liver using bisulfite sequencing in 90 mouse inbred strains, genome-wide expression levels, proteomics, metabolomics, and 68 clinical traits and performed epigenome-wide association studies (EWAS). We found associations with numerous clinical traits including bone density, insulin resistance, expression, and protein and metabolite levels. A large proportion of associations were unique to EWAS and were not identified using GWAS. Methylation levels were regulated by genetics largely in cis, but we also found evidence of trans regulation, and we demonstrate that genetic variation in the methionine synthase reductase gene Mtrr affects methylation of hundreds of CpGs throughout the genome. Our results indicate that natural variation in methylation levels contributes to the etiology of complex clinical traits.

The Genetic Landscape of Hematopoietic Stem Cell Frequency in Mice.

Prior efforts to identify regulators of hematopoietic stem cell physiology have relied mainly on candidate gene approaches with genetically modified mice. Here we used a genome-wide association study (GWAS) strategy with the hybrid mouse diversity panel to identify the genetic determinants of hematopoietic stem/progenitor cell (HSPC) frequency. Among 108 strains, we observed ∼120- to 300-fold variation in three HSPC populations. A GWAS analysis identified several loci that were significantly associated with HSPC frequency, including a locus on chromosome 5 harboring the homeodomain-only protein gene (Hopx). Hopx previously had been implicated in cardiac development but was not known to influence HSPC biology. Analysis of the HSPC pool in Hopx-/- mice demonstrated significantly reduced cell frequencies and impaired engraftment in competitive repopulation assays, thus providing functional validation of this positional candidate gene. These results demonstrate the power of GWAS in mice to identify genetic determinants of the hematopoietic system.

Programmed cell death 1 (PD-1) and its ligand (PD-L1) in common cancers and their correlation with molecular cancer type.

Cancer cells expressing PD-1 ligands (PD-L1/PD-L2) inhibit immune-modulatory T-cell activation facilitating disease progression. Preliminary clinical trials exploring interruption of PD-1/PD-L1 signaling showed benefit in several cancer types. We analyzed the distribution of PD-1-positive tumor-infiltrating lymphocytes (TIL) and cancer cells' expression of PD-L1 in a molecularly profiled cohort of 437 malignancies (380 carcinomas, 33 sarcomas, and 24 melanomas). We showed that the presence of PD-1(+) TILs significantly varied among cancer types (from 0% in extraskeletal myxoid chondrosarcomas to 93% in ovarian cancer), and was generally associated with the increased number of mutations in tumor cells (P = 0.029). Cancer cell expression of PD-L1 varied from absent (in Merkel cell carcinomas) to 100% (in chondro- and liposarcomas), but showed the inverse association with the number of detected mutations (P = 0.004). Both PD-1 and PD-L1 expression were significantly higher in triple-negative breast cancers (TNBC) than in non-TNBC (P < 0.001 and 0.017, respectively). Similarly, MSI-H colon cancers had higher PD-1 and PD-L1 expression than the microsatellite stable tumors (P = 0.002 and 0.02, respectively). TP53-mutated breast cancers had significantly higher PD-1 positivity than those harboring other driver mutations (e.g., PIK3CA; P = 0.002). In non-small cell lung cancer, PD-1/PD-L1 coexpression was identified in 8 cases (19%), which lacked any other targetable alterations (e.g., EGFR, ALK, or ROS1). Our study demonstrated the utility of exploring the expression of two potentially targetable immune checkpoint proteins (PD-1/PD-L1) in a substantial proportion of solid tumors, including some aggressive subtypes that lack other targeted treatment modalities.

Oncocytoma-like renal tumor with transformation toward high-grade oncocytic carcinoma: a unique case with morphologic, immunohistochemical, and genomic characterization.

Renal oncocytoma is a benign tumor with characteristic histologic findings. We describe an oncocytoma-like renal tumor with progression to high-grade oncocytic carcinoma and metastasis. A 74-year-old man with no family history of cancer presented with hematuria. Computed tomography showed an 11 cm heterogeneous multilobulated mass in the right kidney lower pole, enlarged aortocaval lymph nodes, and multiple lung nodules. In the nephrectomy specimen, approximately one third of the renal tumor histologically showed regions classic for benign oncocytoma transitioning to regions of high-grade carcinoma without sharp demarcation. With extensive genomic investigation using single nucleotide polymorphism-based array virtual karyotyping, multiregion sequencing, and expression array analysis, we were able to show a common lineage between the benign oncocytoma and high-grade oncocytic carcinoma regions in the tumor. We were also able to show karyotypic differences underlying this progression. The benign oncocytoma showed no chromosomal aberrations, whereas the high-grade oncocytic carcinoma showed loss of the 17p region housing FLCN (folliculin [Birt-Hogg-Dubé protein]), loss of 8p, and gain of 8q. Gene expression patterns supported dysregulation and activation of phosphoinositide 3-kinase (PI3K)/v-akt murine thymoma viral oncogene homolog (Akt), mitogen-activated protein kinase (MAPK)/extracellular-signal-regulated kinase (ERK), and mechanistic target of rapamycin (serine/threonine kinase) (mTOR) pathways in the high-grade oncocytic carcinoma regions. This was partly attributable to FLCN underexpression but further accentuated by overexpression of numerous genes on 8q. In the high-grade oncocytic carcinoma region, vascular endothelial growth factor A along with metalloproteinases matrix metallopeptidase 9 and matrix metallopeptidase 12 were overexpressed, facilitating angiogenesis and invasiveness. Genetic molecular testing provided evidence for the development of an aggressive oncocytic carcinoma from an oncocytoma, leading to aggressive targeted treatment but eventual death 39 months after the diagnosis.

Genetic regulation of mouse liver metabolite levels.

We profiled and analyzed 283 metabolites representing eight major classes of molecules including Lipids, Carbohydrates, Amino Acids, Peptides, Xenobiotics, Vitamins and Cofactors, Energy Metabolism, and Nucleotides in mouse liver of 104 inbred and recombinant inbred strains. We find that metabolites exhibit a wide range of variation, as has been previously observed with metabolites in blood serum. Using genome-wide association analysis, we mapped 40% of the quantified metabolites to at least one locus in the genome and for 75% of the loci mapped we identified at least one candidate gene by local expression QTL analysis of the transcripts. Moreover, we validated 2 of 3 of the significant loci examined by adenoviral overexpression of the genes in mice. In our GWAS results, we find that at significant loci the peak markers explained on average between 20 and 40% of variation in the metabolites. Moreover, 39% of loci found to be regulating liver metabolites in mice were also found in human GWAS results for serum metabolites, providing support for similarity in genetic regulation of metabolites between mice and human. We also integrated the metabolomic data with transcriptomic and clinical phenotypic data to evaluate the extent of co-variation across various biological scales.

Extraskeletal myxoid chondrosarcoma of the vulva with PLAG1 gene activation: molecular genetic characterization of 2 cases.

Extraskeletal myxoid chondrosarcoma (EMC) is a rare mesenchymal neoplasm, rarely reported in the genitourinary tract with only 5 cases reported in the vulva. We investigated 2 cases of vulvar sarcomas whose morphologic appearance and immunohistochemical profiles were consistent with EMC using fluorescence in situ hybridization (FISH), reverse-transcription polymerase chain reaction, and a whole genome expression array. FISH and reverse-transcription polymerase chain reaction assays showed no EWSR1 and NR4A3 loci rearrangements. Microarray-based analysis also revealed no changes in NR4A3 and EWSR1 gene transcription levels. Microarray data showed a significant downregulation of the muscle-related genes (eg, myosin heavy chain family, actins, myoglobin, desmin, creatine kinase, troponins) and cytokeratins (KRT6A, 6B, 13, 14, and 78), upregulation of several neuron-specific genes [neural cell adhesion molecule 1 (NCAM-1/CD56), neurofilament (NEFH)], along with some well-characterized tumor biomarkers [carbonic anhydrase IX (CA-9), topoisomerase IIα (TOP2A), matrix metalloproteinases (MMP-7, MMP-9), CDKN2 gene (p16-INK4a), checkpoint homolog 2 (CHEK2)]. Notably, both tumors showed upregulation of the pleomorphic adenoma gene 1 (PLAG1), and in 1 case PLAG1 gene rearrangement was detected by break-apart FISH. Some vulvar tumors with morphologic and immunohistochemical characteristics of EMC may represent a molecular genetic entity separate from EMCs arising in other locations. PLAG1 gene activation appears to be involved in the development of these neoplasms.

Analysis of allele-specific expression in mouse liver by RNA-Seq: a comparison with Cis-eQTL identified using genetic linkage.

We report an analysis of allele-specific expression (ASE) and parent-of-origin expression in adult mouse liver using next generation sequencing (RNA-Seq) of reciprocal crosses of heterozygous F1 mice from the parental strains C57BL/6J and DBA/2J. We found a 60% overlap between genes exhibiting ASE and putative cis-acting expression quantitative trait loci (cis-eQTL) identified in an intercross between the same strains. We discuss the various biological and technical factors that contribute to the differences. We also identify genes exhibiting parental imprinting and complex expression patterns. Our study demonstrates the importance of biological replicates to limit the number of false positives with RNA-Seq data.

The systems genetics resource: a web application to mine global data for complex disease traits.

The Systems Genetics Resource (SGR) (http://systems.genetics.ucla.edu) is a new open-access web application and database that contains genotypes and clinical and intermediate phenotypes from both human and mouse studies. The mouse data include studies using crosses between specific inbred strains and studies using the Hybrid Mouse Diversity Panel. SGR is designed to assist researchers studying genes and pathways contributing to complex disease traits, including obesity, diabetes, atherosclerosis, heart failure, osteoporosis, and lipoprotein metabolism. Over the next few years, we hope to add data relevant to deafness, addiction, hepatic steatosis, toxin responses, and vascular injury. The intermediate phenotypes include expression array data for a variety of tissues and cultured cells, metabolite levels, and protein levels. Pre-computed tables of genetic loci controlling intermediate and clinical phenotypes, as well as phenotype correlations, are accessed via a user-friendly web interface. The web site includes detailed protocols for all of the studies. Data from published studies are freely available; unpublished studies have restricted access during their embargo period.

Limited RNA editing in exons of mouse liver and adipose.

Several studies have investigated RNA-DNA differences (RDD), presumably due to RNA editing, with conflicting results. We report a rigorous analysis of RDD in exonic regions in mice, taking into account critical biases in RNA-Seq analysis. Using deep-sequenced F1 reciprocal inbred mice, we mapped 40 million RNA-Seq reads per liver sample and 180 million reads per adipose sample. We found 7300 apparent hepatic RDDs using a multiple-site mapping procedure, compared with 293 RDD found using a unique-site mapping procedure. After filtering for repeat sequence, splice junction proximity, undirectional strand, and extremity read bias, 63 RDD remained. In adipose tissue unique-site mapping identified 1667 RDD, and after applying the same four filters, 188 RDDs remained. In both tissues, the filtering procedure increased the proportion of canonical (A-to-I and C-to-U) editing events. The genomic DNA of 12 RDD sites among the potential 63 hepatic RDD was tested by Sanger sequencing, three of which proved to be due to unreferenced SNPs. We validated seven liver RDD with Sequenom technology, including two noncanonical, Gm5424 C-to-I(G) and Pisd I(G)-to-A RDD. Differences in diet, sex, or genetic background had very modest effects on RDD occurrence. Only a small number of apparent RDD sites overlapped between liver and adipose, indicating a high degree of tissue specificity. Our findings underscore the importance of properly filtering for bias in RNA-Seq investigations, including the necessity of confirming the DNA sequence to eliminate unreferenced SNPs. Based on our results, we conclude that RNA editing is likely limited to hundreds of events in exonic RNA in liver and adipose.

Juxtaglomerular cell tumor: a morphological, immunohistochemical and genetic study of six cases.

Juxtaglomerular cell tumors (JGCTs) are rare tumors characterized by renin synthesis, hyperaldosteronism and hypertension. A curious immunohistochemical overlap between JGCT and gastrointestinal stromal tumor (GIST) including the expression of vimentin, CD34, CD117, α-smooth muscle actin was previously reported, prompting us to further investigate JGCT and its phenotypic and molecular genetic characteristics. Virtual karyotyping showed gain of chromosomes 3, 4, 10, 13, 17 and 18 in one JGCT, and fluorescence in situ hybridization (FISH) study confirmed this multiple gain pattern. Additionally, loss of chromosome 9 was observed in four of six cases analyzed with FISH. A whole genome expression analysis revealed 415 up-regulated (including renin, and CD117) and 325 down-regulated genes between the 2 cases. The study confirmed earlier reports on the gain of chromosomes 4 and 10, and provided further evidence of up-regulation of the genes located on these 2 chromosomes. For the first time our study indicated the importance of the loss of chromosome 9 and loss of expression of several tumor suppressor genes located on this chromosome as possible pathogenetic events important in development of JGCT.

Unraveling inflammatory responses using systems genetics and gene-environment interactions in macrophages.

Many common diseases have an important inflammatory component mediated in part by macrophages. Here we used a systems genetics strategy to examine the role of common genetic variation in macrophage responses to inflammatory stimuli. We examined genome-wide transcript levels in macrophages from 92 strains of the Hybrid Mouse Diversity Panel. We exposed macrophages to control media, bacterial lipopolysaccharide (LPS), or oxidized phospholipids. We performed association mapping under each condition and identified several thousand expression quantitative trait loci (eQTL), gene-by-environment interactions, and eQTL "hot spots" that specifically control LPS responses. We used siRNA knockdown of candidate genes to validate an eQTL hot spot in chromosome 8 and identified the gene 2310061C15Rik as a regulator of inflammatory responses in macrophages. We have created a public database where the data presented here can be used as a resource for understanding many common inflammatory traits that are modeled in the mouse and for the dissection of regulatory relationships between genes.

Hybrid mouse diversity panel: a panel of inbred mouse strains suitable for analysis of complex genetic traits.

We have developed an association-based approach using classical inbred strains of mice in which we correct for population structure, which is very extensive in mice, using an efficient mixed-model algorithm. Our approach includes inbred parental strains as well as recombinant inbred strains in order to capture loci with effect sizes typical of complex traits in mice (in the range of 5% of total trait variance). Over the last few years, we have typed the hybrid mouse diversity panel (HMDP) strains for a variety of clinical traits as well as intermediate phenotypes and have shown that the HMDP has sufficient power to map genes for highly complex traits with resolution that is in most cases less than a megabase. In this essay, we review our experience with the HMDP, describe various ongoing projects, and discuss how the HMDP may fit into the larger picture of common diseases and different approaches.

Genetic locus on mouse chromosome 7 controls elevated heart rate.

Elevated heart rate (HR) is a risk factor for cardiovascular diseases. The goal of the study was to map HR trait in mice using quantitative trait locus (QTL) analysis followed by genome-wide association (GWA) analysis. The first approach provides mapping power and the second increases genome resolution. QTL analyses were performed in a C3HeB×SJL backcross. HR and systolic blood pressure (SBP) were measured by the tail-cuff plethysmography. HR was ∼80 beats/min higher in SJL compared with C3HeB. There was a wide distribution of the HR (536-763 beats/min) in N2 mice. We discovered a highly significant QTL (logarithm of odds = 6.7, P < 0.001) on chromosome 7 (41 cM) for HR in the C3HeB×SJL backcross. In the Hybrid Mouse Diversity Panel (58 strains, n = 5-6/strain) we found that HR (beats/min) ranged from 546 ± 12 in C58/J to 717 ± 7 in MA/MyJ mice. SBP (mmHg) ranged from 99 ± 6 in strain I/LnJ to 151 ± 4 in strain BXA4/PgnJ. GWA analyses were done using the HMDP, which revealed a locus (64.2-65.1 Mb) on chromosome 7 that colocalized with the QTL for elevated HR found in the C3HeB×SJL backcross. The peak association was observed for 17 SNPs that are localized within three GABA(A) receptor genes. In summary, we used a combined genetic approach to fine map a novel elevated HR locus on mouse chromosome 7.

Quantitative trait Loci association mapping by imputation of strain origins in multifounder crosses.

Although mapping quantitative traits in inbred strains is simpler than mapping the analogous traits in humans, classical inbred crosses suffer from reduced genetic diversity compared to experimental designs involving outbred animal populations. Multiple crosses, for example the Complex Trait Consortium's eight-way cross, circumvent these difficulties. However, complex mating schemes and systematic inbreeding raise substantial computational difficulties. Here we present a method for locally imputing the strain origins of each genotyped animal along its genome. Imputed origins then serve as mean effects in a multivariate Gaussian model for testing association between trait levels and local genomic variation. Imputation is a combinatorial process that assigns the maternal and paternal strain origin of each animal on the basis of observed genotypes and prior pedigree information. Without smoothing, imputation is likely to be ill-defined or jump erratically from one strain to another as an animal's genome is traversed. In practice, one expects to see long stretches where strain origins are invariant. Smoothing can be achieved by penalizing strain changes from one marker to the next. A dynamic programming algorithm then solves the strain imputation process in one quick pass through the genome of an animal. Imputation accuracy exceeds 99% in practical examples and leads to high-resolution mapping in simulated and real data. The previous fastest quantitative trait loci (QTL) mapping software for dense genome scans reduced compute times to hours. Our implementation further reduces compute times from hours to minutes with no loss in statistical power. Indeed, power is enhanced for full pedigree data.